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I. INTRODUCTION 

Transverse densities have become an essential tool in the analysis of current matrix elements (vector, axial) and 
the description of the spatial structure of hadrons M- They are defined as the two-dimensional Fourier transforms 
of the invariant form factors and describe the distribution of charge and current in the hadron in transverse space 
at fixed light-front time x + = t + z. They are frame-independent (boost-invariant) and provide an objective spatial 
representation of the hadron as a relativistic system. In the context of QCD the transverse densities correspond to a 
projection of the generalized parton distributions (GPDs) describing the transverse spatial distribution of quarks and 
antiquarks ii; as such they connect the information gained from low-energy elastic scattering with the partonic 
content probed by high-momentum-transfer processes in high-energy scattering. In composite models of nucleon 
structure the transverse densities can be expressed as proper densities of the light-front wave functions of the system 
and are therefore a natural ground for phenomenological analysis. Considerable efforts have been devoted to extracting 
the transverse charge and magnetization densities in the nucleon from the available electromagnetic (Dirac, Pauli) 
form factor data and studying their properties d, @, 01 ; see Ref. [I[ for a review. 

Of particular interest are the densities in the nucleon’s chiral periphery. At transverse distances b = 
where the pion mass is regarded as parametrically small compared to the typical inverse hadronic size, the densities 
are governed by the universal dynamics resulting from the spontaneous breaking of chiral symmetry and can be 
computed from first principles using the methods of chiral effective field theory (EFT). The isovector charge and 
magnetization densities, p\{b) and P 2 (b), arise from chiral processes in which the current couples to the nucleon 
through exchange of a two-pion system with momenta 0{M„) relative to the nucleon. A detailed investigation of 
the properties of the peripheral transverse densities in leading-order (LO) relativistic chiral EFT was performed in 
Ref. [sj. The densities decay exponentially with a range given by the mass of the exchanged system, 2 M n (“Yukawa 
tail”); their overall strength and the underlying power-like behavior in b are determined by coupling of the two-pion 
exchange to the nucleon and exhibit a rich structure. It was found that the transverse charge density p\ (b) and the 
modified magnetization density p%(b) = (d/db)pX (6)/(2Mjy) are of the same order in the chiral expansion, obey an 
approximate inequality pX (b) < Pi (b), and are numerically very close at the distances of interest, b > 1 M ~ 1 . These 
findings represent model independent features of the nucleon’s chiral periphery and call for a simple explanation. 

In Ref. [§] the peripheral densities were calculated in a dispersive representation, where they are expressed as 
integrals of the imaginary parts (or spectral functions) of the invariant form factors along the cut in timelike region at 
t > 4 M%. This formulation makes it possible to use the well-known chiral EFT results for the invariant form factors 
and their spectral functions for the calculation of the transverse densities. While it allows one to derive all properties 
of interest, it does not provide a mechanical picture of the chiral processes as pions “moving about” the nucleon in 
space and time. Such a picture could be obtained in a time-ordered representation of chiral EFT, where one works 
with the concepts of instantaneous configurations, time evolution, and the wave function of the system. Since the 
transverse densities are defined at fixed light-front time x + it is natural to adopt light-front quantization nnm and 
follow the evolution of the relevant chiral processes in light-front time. This representation might explain our earlier 
findings and provide new insight into the structure of the peripheral transverse densities. 

Studying the space-time evolution of chiral dynamics in light-front quantization is interesting also for methodolog¬ 
ical reasons, unrelated to the specific questions posed by transverse densities. The typical momentum of soft pions in 
the nucleon rest frame is k = 0{M. K ) [the velocity is v = 0(1)], and the typical energy of configurations is E = 0(M n ). 
Chiral dynamics thus represents an essentially relativistic system, in which pions “appear” and “disappear” through 
quantum fluctuations and the number of particles is generally not conserved. In equal-time quantization the particle 
number observed at an instant changes under Lorentz boost, so that the wave function is essentially frame-dependent 
and no meaningful particle-based description of the theory can be constructed. In light-front quantization the particle 
number is invariant under boosts, the wave function is frame-independent, and a natural particle-based description 
is obtained EH- It represents the only known formulation that permits a consistent first-quantized particle-based 
description of chiral processes. Such a representation could significantly advance our understanding of chiral dynamics. 

In this article we study the nucleon’s transverse charge and magnetization densities in the chiral periphery in a 
first-quantized particle-based representation of chiral dynamics based on light-front quantization. The LO chiral 
EFT results for the peripheral densities are expressed in time-ordered form, as the result of a transition of the bare 
nucleon to a virtual irN state mediated by the chiral EFT interactions. The densities appear as overlap integrals of the 
perturbative light-front wave functions describing the N —> ttN transition, which are calculable directly from the chiral 
Lagrangian. The new representation offers new insight into the structure of peripheral densities and reveals several 
interesting properties. First, it explains in simple terms the parametric order of the peripheral charge and modified 
magnetization densities, p\{b) and p 2 (b), in the chiral expansion. Second, it proves the inequality \pX (b)\ < p\{b ), 
which had been observed numerically in the earlier study using the invariant formulation [Sj. It also explains why 
the inequality is approximately saturated, pX (b) ss p\{b), and shows that this is related to the essentially relativistic 
character of chiral dynamics. Third, the wave function overlap representation exposes the role of the pion’s orbital 


4 


angular momentum in the peripheral transverse densities. A particularly simple picture is obtained with transversely 
polarized nucleon states, where only a single pion orbital with L = 1 accounts for both densities, and the relation 
between p \ (&) and P\ (b) is explained by the “left-right” asymmetry induced by the orbital motion of the pion in the 
preferred longitudinal direction. We emphasize that the first-quantized representation developed in the present work 
is equivalent to the invariant LO chiral EFT expressions used in earlier studies, and that the new insights derived 
from it reflect general properties of the nucleon’s chiral periphery. A summary of our results has been presented in 
Ref. jjj]. 

In the present work we derive the light-front representation of peripheral chiral processes by rewriting the result 
obtained in Lorentz-invariant chiral EFT. This approach has the advantage that it guarantees equivalence to the 
well-tested invariant formulation and avoids the use of light-front specific techniques. The nN light-front wave 
function of the nucleon in chiral EFT is defined in terms of the vertex function provided by the chiral Lagrangian, 
and the wave function overlap representation of the current matrix element is obtained naturally from the reduction 
of the Feynman integrals. The connection with the conventional light-front time-ordered Hamiltonian approach El 
is explained in Appendix [A] Some formal aspects of chiral EFT in the time-ordered formulation were studied in 
Refs. [TM3. Our derivation also confirms the presence of an instantaneous term (or zero mode contribution) in the 
light-front representation of the chiral component of the transverse charge density 0®. This term has a simple 
physical interpretation as describing the contribution of large-mass (non-chiral) intermediate states in time-ordered 
chiral processes and is shown to be numerically small. 

The light-front representation of chiral dynamics described here can be applied also to the peripheral transverse 
densities of other local operators, e.g. the matter and momentum densities associated with the energy-momentum 
tensor, and to the peripheral GPDs probed in high-energy scattering processes. We derive the wave function overlap 
representation of the light-front “plus” momentum density of peripheral pions in the nucleon, which determine the 
chiral component of the nucleon’s parton densities at transverse distances b = 0(M~ 1 ). This establishes a formal 
connection between our chiral EFT results for the peripheral transverse densities and the nucleon’s quark/antiquark 
content. The light-front momentum density of peripheral pions could in principle also be probed directly in peripheral 
high-energy scattering processes. 

The plan of this article is as follows. In Sec. m we review the transverse density representation of the current 
matrix element, the peripheral chiral contributions in the invariant formulation [$], and derive the light-front overlap 
representation of the current matrix element. In Sec. lIYI we investigate the properties of the peripheral ttN light-front 
wave function, including the choice of nucleon spin states and the coordinate representation. In Sec.|V] we express 
the transverse densities p\ and p 2 as overlap integrals of the coordinate-space nN light—front wave functions and 
study their properties. Using longitudinal nucleon spin states we discuss the parametric order of the densities, derive 
the inequality between them, and evaluate them numerically. We also present the expressions for transverse nucleon 
spin states and show that they correspond to a simple mechanical picture of a pion with L = 1 orbiting around the 
nucleon in the rest frame. This picture concisely summarizes the dynamical content of the LO chiral EFT contribution 
and represents the main result of this work. We also compute the instantaneous (contact term) contribution to the 
densities and show that it is numerically small. In Sec. El we connect the transverse densities with the peripheral 
parton content of the nucleon in QCD. We derive the wave function overlap representation of the pion plus momentum 
distribution (“pion GPD”) in chiral EFT, which determines the nucleon’s peripheral parton densities ®, ®|, and 
show that the transverse charge density is recovered by integrating the peripheral quark/antiquark densities over the 
parton momentum fraction x. 


II. TRANSVERSE DENSITIES 


The transition matrix element of the electromagnetic current between nucleon states is parametrized in terms of 
two invariant form factors (we follow the notation and conventions of Ref. Q) 


{N{p 2l a 2 )\J^{x)\N{jp ll a 1 )) = u 2 




u i e 


„iA.x 


( 2 . 1 ) 


where pi, 2 are the nucleon 4-momenta, (71,2 the spin quantum numbers, u\ = u(.Pi, cr i) etc- the nucleon bispinors, 
normalized to uiui = u 2 u 2 = 2A/,y, and a t [ 7 ^, 7 "]- The 4-momentum transfer is defined as 

A = p 2 -pi, (2.2) 


and the dependence of the matrix element on the space-time point x where the current is measured is dictated by 
translational invariance. The Lorentz-invariant momentum transfer is 


= A 2 = {p 2 -pi) 2 , 


t 


(2.3) 
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FIG. 1. Interpretation of the transverse densities of the electromagnetic current in a nucleon state with spin quantized in 
the transverse ^-direction, Eq. (ED) . pi(b) describes the spin-independent (or left-right symmetric) part of the plus current 
density; cos<j>p 2 (b) describes the spin-dependent (or left-right asymmetric) part. 


with t < 0 in the physical region for electromagnetic scattering. The Dirac and Pauli form factors, Fj and F 2 , are 
invariant functions of t and can be discussed independently of any reference frame. 

In the context of the light-front description of nucleon structure one naturally considers the form factors in a frame 
where the momentum transfer vector lies in the transverse ( x-y ) plane, 


= (A 0 , A X ,A V , A z ) = (0, A t , 0), A T = (A x ,A y ), 
and represents them as Fourier transforms of certain spatial densities 31 


t = - A 




F\^{t — ~^r) — 


d 2 b e iAT ' b 


Pi, 2 ( 6 ), 


(2.4) 


(2.5) 


where b = (b x ,b v ) is a transverse coordinate variable and b = \b\. The formal properties of the transverse densities 
Pi, 2 (b) and their physical interpretation have been discussed extensively in the literature 31 - and are summarized 
in Ref. [§]. They describe the transverse spatial distribution of the light-front plus component of the current, J + = 
J° ± J z , in the nucleon at fixed light-front time x + . Specifically, in a state where the nucleon is localized in transverse 
space at the origin, and polarized in the y-direction, the matrix element of the current J + at light-front time x + = 0 
and light-front coordinates x~ = 0 and x T = b is given by 

(J + (6))iocaiized = (■•■) [pi(b) + (2S V ) C0S(f> p 2 (b)} , (2.6) 


k(4) s ¥b 


P2(b) 


2Mjv 


(2.7) 


where (...) hides a trivial factor reflecting the normalization of states (see Ref. [8i] for details); cos 4> = b x /b is the 
cosine of the azimuthal angle, and S v = ±1/2 the spin projection in the y-direction in the nucleon rest frame (see 
Fig. HD- The function p\(b) describes the spin-independent part of the current; the function cos (f>p 2 (b) describes the 
spin-dependent part of the current in a transversely polarized nucleon. 

The spin-dependent part of the current Eq. (12.611 changes sign between negative and positive values of b x , or “left” 
and “right” positions when looking at the nucleon along the negative z-direction (down from z = ±oo; see Fig. [TJ. 
The density p 2 (b) can thus be interpreted as the left-right asymmetry of the J + current in a nucleon polarized in the 
positive y direction. This interpretation has a natural connection with composite models of nucleon structure, where 
polarization in the y-direction generally induces a convection motion of the constituents around the y-axis. As a 
result, the observer sees the charged constituents on the left side as “blue-shifted” (larger plus momentum), and those 
on the right ride as “red-shifted” (smaller plus momentum), compared to the unpolarized case Q. We shall refer 
to this interpretation in our discussion of the peripheral chiral component in Sec. IVDl We note that Eq. (12.61) and 
its interpretation can be generalized to the case of arbitrary nucleon polarization states in the rest frame, including 
non-diagonal transitions; see Ref. JJ for details. 
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The electromagnetic current matrix element and the transverse densities have two isospin components. In the 
following we are concerned with the isovector component 

(iV|J|A) y = \ [{p\J\p) - (n\J\n )], p y 2 = |(Pi, 2 “ Pi, 2 )- ( 2 -8) 

The isoscalar component is defined by the same expression with the + sign. 


III. CHIRAL DYNAMICS IN CURRENT MATRIX ELEMENT 


A. Peripheral chiral processes 


At distances b = 0(M V 1 ) the transverse densities are governed by universal chiral dynamics and can be calculated 
using methods of chiral EFT Id- The isovector charge and magnetization densities in this region arise from chiral 
processes in which the current couples to the nucleon through exchange of a two-pion system in the t-channel. At 
LO these are the processes described by the Feynman diagrams of Fig. [5Jr, where the vertices denote the pion nucleon 
couplings of the LO relativistic chiral Lagrangian 0]. They produce densities of the form 

p£ 2 (&) ~ Pi, 2 (M„,M N -b) exp(-2M 7r 6), (3.1) 


where the the exponential decay is determined by the minimal mass of the exchanged system, 2Af 7r ; the pre-exponential 
factors Py 2 is determined by the coupling of the exchanged system to the nucleon and exhibits a rich structure due 
to its dependence on the two scales, M n and Mjy. Diagrams in which the current couples directly to the nucleon, or 
to a pion-nucleon vertex, produce contributions to the densities with range 0(M^ 1 ) Lor terms oc 6^ 2 \b), and do not 
need to be considered in the calculation of the densities at b = 0(AL“ 1 )|J In Ref. [§[ the densities at b = 0(M~ 1 ) 
resulting from the diagrams of Fig. [3^, were computed in a dispersive representation, where the densities are expressed 
as integrals of the imaginary parts of the form factors on the cut at t > 4Mi;, and the integral extends over the 
parametric region t — AM 2 = 0{M 2 ). 

Now we want to represent the chiral dynamics generating the peripheral densities as actual processes evolving 
in light-front time x + , and to express the densities in terms of the light-front wave functions of the chiral nN 
system. This could be done by solving the dynamical problem of chiral EFT directly using light-front time-ordered 
perturbation theory 0- A more convenient approach is to take the known chiral EFT result in the relativistically 
invariant formulation and rewrite it such that it corresponds to the overlap of light-front wave functions. This 
approach maintains the connection with the invariant formulation and naturally generates also the instantaneous 
terms (zero modes) that require special considerations in the time-ordered approach. In the LO approximation the 
invariant chiral EFT result for the isovector nucleon current matrix element Eq. ED at the position x = 0 is 0-0] 
(the specific form here was derived in Ref. Q) 


{N 2 \ ^(0) |Vi) y 


d 4 k i 

(27r) 4 {k 2 - Ml + i0){k% - Ml + i0) 


9 a ^2*275^ + M N )kij 5 ui 
Fl l 2 - M 2 n + *0 


+ -^2 U 2*«l 


where 


+ (diagrams without 7T7 t cut), 


(3.2) 


* 1,2 = kT A/2 (3.3) 

are the 4-momenta of the pions coupling to the vector current, and the average 4-momentum k was chosen as 
integration variable. The first term in the brackets results from the triangle diagram in Fig. [5K and is proportional 
to the squared nNN coupling g\/F 2 \ it involves the intermediate nucleon propagator with 4-momentum 

l=p 1 -ki = p 2 - k 2 (3.4) 

(we use the notation * lj2 = ^ — ^7/J- The second term results from the contact diagram in Fig. and is 

proportional to the mrNN contact coupling in the chiral Lagrangian, 1 /F 2 . As explained above, Eq. (13.21) shows only 
the contribution from the rnr cut diagrams that contribute to te peripheral density. 


1 The diagrams in which the current couples directly to the nucleon or to a pion-nucleon vertex renormalize the charge in the center of 
the nucleon, to compensate for the peripheral charge density produced by the two-pion exchange diagrams and ensure overall charge 
conservation. 
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^2 ^ ^ + k 2 l'"\ .;'"t + (no 7I7T cut) 

* \ \ * 


■ ■ * - ■ ■ < < 

P 2 l Pi Pi Pi 

(a) 


(b) 


& 


,< 8 K 


+ 


t t 

8a Ys / 


t 

(1 ~8a) IFn 


+ (no 7171 cut) 


(C) 


j A T - P2j~P 


IT 


>’Pt 


( 8 ) 


& T + p 


/ 


2T / 

/ 


ypl 

s ^ ^T+ZDt 

\ 


P\ P 2T (1-V)P + , “*T 


P\ P IT 


FIG. 2. (a) Feynman diagrams of LO chiral EFT processes contributing to the the two-pion cut of the isovector nucleon form 
factors and the peripheral transverse densities, Eq. (GH). The squares denote the vertices of the relativistic chiral Lagrangian 
(axial vector nNN coupling, nnNN contact coupling), (b) Equivalent representation Eqs. (13.711 (13.91 . The circles denote 
the pseudoscalar nNN coupling and the effective contact coupling oc (1 — g\). (c) Light-front representation of the triangle 
diagram. The original nucleon makes a transition to a pion-nucleon intermediate state that couples to the current. The plus 
and transverse momenta are indicated. 


The first term in Eq. (lira contains a piece in which the pole of the nucleon propagator cancels, and which is of 
the same form as the second term. For deriving a wave function overlap representation it is important to extract 
this piece and combine it with the second term. Expressing the pion momenta as k 1,2 = Pi ,2 — l, cf. Eq. m, using 
the anticommutation relations between the gamma matrices, and the Dirac equation for the external nucleon spinors, 
Piiii = Mnu 1 and U 2 P 2 = W 2 -Mat, we rewrite the bilinear form in the numerator of the first term in Eq. m as 

^2^275 (l + M N )ki 75U1 = -4M^U275(i + Mat)75Wi - (l 2 - Mlf)[u 2 kui + (terms even in A:)]. ( 3 . 5 ) 

In the second term the factor l 2 — M^ cancels the pole of the nucleon propagator. Moreover, the terms even under 
k —t—k in the parentheses integrate to zero because after canceling the nucleon pole the remaining integrand in 
Eq. (13.21) is odd under k -A —k. Altogether, we can thus replace the terms in the bracket in Eq. (13.21) by 


4 m n9a u 2ls{l + M n )^Ui 
FI P-M 2 n + zO 


1 ~9a 

n 


U2ku\ 


(3.6) 


The peripheral chiral EFT contribution to the current matrix element can therefore be represented as the sum of two 
terms (see Fig. m 

(N 2 \ J»mNl) V = <--->Xterm + (• • ■) contact ; (3-7) 


<iv 2 | J^(0) \N X ) 

int 


{N 2 \ ^(0)1^1) 


am n9a [ 

n 

1 -9a 


(-*) kV U 2 Js{l + Mn)15Ui 


(2?r) 4 (jfe? - Ml + i0)(k 2 - Ml + zO)(Z 2 - M 2 + *0) ’ 
d A k ik 11 v,2ku\ 


contact — 


Fl J (2tt) 4 (k 2 -Ml+iO){k 2 -Ml+iO)' 


(3.8) 

(3.9) 


The first term, Eq. (13.81) , contains the intermediate nucleon propagator and is identical in form to the pion loop graph 
with the usual pseudoscalar nNN vertex with effective coupling g^NN = Mj^qa/F n . This will allow us to derive a 
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wave function overlap representation for this term with the pseudoscalar vertex, which is free of the ambiguities of 
the momentum-dependent axial vector coupling. The second term, Eq. m , represents an effective contact term, 
combining the explicit nirNN 4-point vertex in the chiral Lagrangian with the “non-propagating” part of the triangle 
diagram. The appearance of the combination 1 —g\ indicates that this term expresses internal structure of the nucleon 
(for a pointlike Dirac fermion g J 4 = 1) and that its contribution is numerically small; cf. the discussion in Sec. IV FI 
The two terms in the current matrix element thus have distinct physical meaning and will be discussed separately in 
the following. We emphasize that the decomposition Eqs. (t3TD m is obtained by identical rewriting of the original 
Feynman integrals and does not involve additional approximations. 


B. Overlap representation 

The intermediate-nucleon term of the current matrix element can be represented as an overlap integral of light-front 
wave functions. We derive this representation through a suitable three-dimensional reduction of the Feynman integral 
Eq. (13.811 . To this end we go to a class of frames where the momentum transfer has only transverse components, cf. 
Eq. (12.41) . such that (see Fig. [2t) 


Pi = P2 = 


P 


P2T — Pit = A 


+ Pit „- _ 1V± N t t> 2 T /, in \ 

T, Pi — T ) P 2 ^ + ■ ( 3 .Ill) 

p + p + 

The plus component p + > 0 is a free parameter, whose choice selects a particular frame in a class of frames related 
by longitudinal boosts. Likewise, the overall transverse momentum remains unspecified; only the difference P 2 X — Pit 
is required to be equal to the momentum transfer A x- We introduce light-front components of the loop momentum 
k± = k° ± k z and kx = ( k x , k v ), 

f d A k = \ [ dk+ [ dk~ [ d 2 k T , (3.11) 


m 2 n - 


p\t 


and express k + in terms of the boost-invariant pion momentum fraction y, 

k + = yp + (3. 12 ) 

The integrand of Eq. (13.81) has simple poles in fc - , at the values determined by the mass shell conditions for the pion 
and nucleon 4-momenta. The two poles of the pion propagators and the one of the nucleon propagator lie on opposite 
sides of the real axis if 0 < y < 1. The integral over k~ can thus be taken by closing the contour around the nucleon 
pole. At the nucleon pole the pion virtualities take the values 


k 2 - Ml = - 


(k T + VPit ) 2 V 2 M 


-Ml < 0 , 


(3.13) 


k 2 - Ml = - 


( k T + yP 2 T) 2 y 2 M 


N 


-Ml < 0 , 


(3.14) 


where 


y = 1 -y- 


(3.15) 


These virtualities can be related to the invariant mass differences between states in the light-front time-ordered 
formulation, in which the external nucleon makes a transition to an intermediate nN state and back (see Fig. [2t). 
The invariant mass difference for the transition from the initial nucleon state with plus momentum p + and transverse 
momentum Pit to a pion with yp + and kx + Pit and a nucleon with yp + and — kx is given by 

{kx + Pit ) 2 + Ml 


AM (y, kx, Pit) — 


k t 


M N 


y 


y 


— M^f — p 1T 


(3.16) 


{k T + ypix) 2 + Ml : {k T + ypix) 2 + 


y 


+ 


y 


-m 2 n . 


(3.17) 


The invariant mass difference for the transition from the final nucleon state with p + and j> 2 t to a pion with yp + and 
kx + P2T and a nucleon with yp + and — kx is given by the same expressions with pix —> P2T- It is easy to see that 


k 2 i ~ Ml 

y 

k 2 - Ml 


= AM 2 {y,k T ,pix), 


= AM 2 {y,k T ,P2T)- 


(3.18) 


y 


(3.19) 
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Equations (13.181) and (13.191) allow us to interpret the pion propagators in the Feynman integral as invariant mass 
denominators. [The origin of the expression Eq. (13.161) and its connection with the “energy denominator” in light- 
front time-ordered perturbation theory are explained in Appendix(Aj this information is not needed for the calculations 
performed here but important for general understanding.] 

Further, at the pole of the nucleon propagator the numerator of Eq. (13.81) can be factorized. Since at the pole the 
4-momentum l is on the mass shell, the matrix l + M.y coincides with the projector on physical nucleon spin states 
and can be represented as 

(Mat + 0 on—shell — E u(l,a)u(l,cr), (3.20) 

(7—±1/2 

where u(l, a) is a set of nucleon 4-spinors; the choice of polarization states will be specified below. We thus can write 
the bilinear form in the numerator as [reverting to the full notation u\ = u(pi,ai) and u 2 = u(p 2 ,<J 2 )] 

- M275(^+ Mjv)7 5 Ui [on—shell = ^ u(p 2 , CT 2 )*7 5 u(^ Cr) u(l , a)*7 5 u(pi, (Ji). (3-21) 

<t=±1/2 

Here it is understood that the on-sliell 4-momentum l is expressed in terms of the remaining integration variables y 
and kr, 


1+ = yp + i 


\k T \ 2 + 

yp+ 


It = — A it • 


(3.22) 


The bilinear forms appearing on the right-hand side of Eq. (13.211) can be related to the vertex functions for anlV-> nN 
transition with specified on-shell nucleon momenta and spin and its complex conjugate. Defining the pseudoscalar 
vertex function for the transition from the initial nucleon state with momentum p + and pit and spin a± to a nucleon 
with momentum yp + and —fey and spin a as 


r(2/,fc T ,pi T ;cr, CTi) 


g A M N 


u(l,a)i'y 5 u(pi,a 1 ), 


(3.23) 


the vertex for the transition to the final state is given by 


g A M N M (p 2)0 - 2 )^y 5M (f )0 -) = (7)i 75 u(}, 2i (7 2 )]* = r*(j(, A t ,P2t; ff 2 ), 

F 7 T F TT 


(3.24) 


and multiplying Eq. (13.211) by the squared coupling constant we obtain 


g\M% 

n 


U2"/5{1 + M N )"f 5 Ul\ on—shell 


E r *(2b P 2T ! CT ’ ct 2 )T (y, k T , Pit ; cr,ai). 

<r=± 1/2 


We now define the light- front wave function of the initial state in the process of Fig. \2b as 


^(y, k T ,PiT;cr, a i) 


r(y, fcr. Pit', cr, cri) . 
AM' 2 (y 1 k T ,PiT) 


(3.25) 


(3.26) 


the wave function for the final state is given by the same expression with pit -A P 2 T and a± -A <72 (i.e., it is the same 
function but evaluated at a different argument). It is then straightforward to compute the integral over A: - , and the 
intermediate-nucleon part of the current matrix element Eq. m becomes 

(IVal J + (0) | Ab ) interm = ( N (P + , P2T, CT 2 ) | J+(0) | N ( P + , p 1T , ))L erm 

= J ( 3 - 27 ) 


The original Feynman integral is represented as an overlap integral of the light-front wave functions describing the 
transition from the initial nucleon state N\ to a nN intermediate state and back to the final nucleon state W- 
Equation (13.271) will be our starting point for the analysis of the transverse densities. 

The wave function Eq. (13.261) is defined in terms of the vertex function obtained from the chiral Lagrangian 
and the invariant mass difference of the N -A nN transition. Appendix [A] shows that this object is identical to the 
traditional light-front wave function, defined as the transition matrix element between the initial nucleon state and the 
intermediate pion-nucleon state in light-front time-ordered perturbation theory. Regarding isospin the wave function 
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Eq. (13.26[) is normalized such that it describes the transition p —> 7r° + p, for which the coupling is g^o pp = g^NN = 
gA^N/F n . For the transitions p —> n + +n and n —> n~ +p , which actually contribute to the isovector electromagnetic 
current matrix element, the couplings follow from isospin invariance and are g^+np = 9n~pn = V^g^NN- The isospin 
factor \/2 x \/2 = 2 is included in the prefactor of Eq. (13.271) . 

As explained above, we are interested in the current matrix element only in the (unphysical) region of momentum 
transfers in the vicinity of the two-pion threshold in the t-channel, t — AM% = 0(M 2 ), and consider the wave function 
representation Eq. (13.271) only in this parametric domain. The restriction to this domain will appear naturally when 
going over to the coordinate representation and considering the region of transverse distances b = 0(M~ 1 ). 


IV. CHIRAL LIGHT-FRONT WAVE FUNCTION 


A. Nucleon spin states 


To proceed with the evaluation of the overlap formula Eq. (13.271) we need to specify the nucleon spin states and 
obtain explicit expressions for the vertex functions Eq. (13.231) and (13.241) . Equation (13.271) can be evaluated with any 
choice of nucleon spin states (external and internal); the resulting expressions and their interpretation depend on the 
choice, of course. It is natural to choose the nucleon spin states as light-front helicity states [ll|. For a nucleon state 
with light-front momentum p + and pt the light-front helicity spinors are obtained by subjecting the Dirac spinors 
in the rest frame [p + (RF) = Mm, Pt(RF) = 0] first to a longitudinal boost from Mm to p + , and then to a transverse 
boost from 0 to pr- The spinors thus defined are invariant under longitudinal boosts and transform in a simple 
manner under transverse boosts. An explicit representation of the light-front helicity spinors is Elm 


u(p, cr) = u(p + 1 p T ,cr) 


7 —- [p + 7 + ( m n - IT-PT )7 + ] 

V 2 p 



(4.1) 


where p denotes the on-shell 4-momentum vector ( p 2 = M'^), a = ±1/2, and x(cr) are rest frame 2-spinors for 
polarization in the positive and negative ^-direction, 


X (cr = 1/2) = 



X(cr = -1/2) = 



(4.2) 


The spinors are normalized such that uu = 2 Mm and satisfy the completeness relation Eq. (13.201) . Using these spinors 
to evaluate the pseudoscalar vertex Eq. (13.231) one gete0 

F(y,fc T , pi r ; cr,cri) = [ y M N S z {a, or) ± (k T + VPit) ■ Sricr.ax)] , (4.3) 

FirVV 

and similarly for the vertex with p 2 T and 02 . Here S z and St = (S X ,S V ) are the components of the 3-vector 
characterizing the spin transition matrix element in the rest frame 

S\a,a r) = X t (o-)(|o- l )x(<7 1 ) (■i = x,y,z ), (4.4) 

where a 1 are the Pauli matrices. Use of this compound variable results in a compact representation of the light-front 
spin structure in close correspondence to non-relativistic quantum mechanics. 

The vertex Eq. COD contains two structures with different orbital angular momentum (see Fig. ED- The first term 
on the right-hand side is diagonal in the light-front helicity, because 

S z (<T,<Ti) = a S(cr, a x ) (4.5) 

when the rest frame spinors are eigenspinors of cr 2 , cf. Eq. (14.21) . It describes a transition N —> ttN in which the 
nucleon light-front helicity is preserved and the nN state has orbital angular momentum projection L z = 0. The 
second term is off-diagonal in light-front helicity, because a x and cr v have only off-diagonal elements. It corresponds 
to a transition in which the nucleon helicity is flipped and the nN system has orbital angular momentum projection 
L z = 1. This is immediately obvious from the fact that this term is proportional to the transverse momentum 
kr ± ypiTi which transforms as 2-dimensional vector under rotations around the z-axis. 


2 The sign of the three-dimensional expressions for the vertex function depends on the convention for the matrix 7 ®. We use the 
Bjorken-Drell convention 7 ® = 27 ° 7 1 7 2 7 ®. 
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FIG. 3. (a) Spin structure of the N —» nN light-front wave function in chiral EFT. The light-front helicity corresponds to 
the ^-projection of the nucleon spin in the rest frame. The helicity-conserving component has orbital angular momentum 
projection L z = 0, the helicity-flip component has L z = 1. (b) Wave function overlap representation of the transverse densities 
pY (b) and pY {b) in the light-front helicity representation, Eq. (15.411 . 


B. Transverse rest frame 


Equation (13.271) represents the current matrix element as an overlap integral of the light-front wave functions of the 
initial and final nucleon states with overall transverse momenta Pit and p 2 T- For further analysis is will be convenient 
to express these wave functions in terms of the wave function at zero overall transverse momentum (transverse rest 
frame), such that the overlap integral becomes a quadratic form in a single function. This can be accomplished using 
the transformation properties under transverse boosts. As can be seen from Eqs. (13.261) . (13.171) and (14.31) . the wave 
function at overall transverse momentum pit is related to that at zero transverse momentum by 


= ^(y,k T + ypiT, 0;cr,cri) = Hf(y, k T + ypir] a, d). 


(4.6) 


The expression for the rest frame wave function can be obtained by setting p\r = 0 in Eqs. (13.171) and (14.31) . For 
reference we quote the explicit formulas: 


^(y,k T ;a, cxi) = 


A M 2 (y, k T ) = — 


T (y,k T ;cr, (Ji) 
A M 2 {y 1 k T ) 

Ml 


kly 


T(j/, k T ] ct, <7i) = 


y 

2igAM N 

FttVv 


' M N _ M 2 
- m Ni 


yM N S z (a , (7i) + k T ■ S T (cr, ai) 


(4.7) 

(4.8) 

(4.9) 


where we use fcr to denote the transverse momentum argument. Similar formulas apply to the outgoing wave function 
with transverse momentum p 2 T and spin (72. The current matrix element Eq. (13.271) can therefore equivalently be 
expressed in terms of the rest frame wave function, 


(iV 2 | J+(0) \N,) 

i 


V 


(2 P 4 


dy_ 
yy , 


dlkx 


2 tt J yy J {2 i r) 2 ^ 

We shall use this expression in our theoretical studies in the following. 


^2 k T + yP 2T ’ a > CT 2) ^(y, k T + yPlT', cr, CTl). 


(4.10) 


C. Coordinate representation 

It is instructive to study the rest frame light-front wave function in the transverse coordinate representation. 
The coordinate-space wave function allows us to identify the parametric regime of peripheral distances where chiral 
dynamics is valid and to calculate the transverse densities directly in coordinate space. We define the coordinate-space 
wave function as the transverse Fourier transform of the momentum-space wave function at fixed plus momentum 
fraction y, 


(2tt)2 


IVT’^T 


T(y,fc T ;cr, (Ti). 


®{y,r T ,cr,o- 1 ) 


1.1 1 



























12 


The vector rx is the difference in the transverse positions of the n and N (relative transverse coordinate), such that 
the wave function describes the physical transverse size distribution of the nN system. Because of the spin structure 
of the vertex Eq. S3) the coordinate-space wave function can be expressed in terms of two transverse radial wave 
functions (i.e., scalars with respect to rotations around the 2 -axis), 


$(y,r T ,cr, err) = -2iS z (a,ai) U 0 (y,r T ) + 


2 r T ■ S T {a , <Ti) 


r T 


Ui{y,r T ), 


(4.12) 


where rx = \rx\ is the modulus of the transverse coordinate. Following Sec. IIV Al these are the components with 
orbital angular momentum projection L z = 0 and 1. The Fourier integral is easily calculated by writing the invariant 
mass in the denominator of the momentum-space wave function, Eq. ()4.8I) . in the form 


k 2 

AM 2 = T 


■M 2 


yy 


M t = M T (y) = ^JyM 2 +y 2 M 2 f , 

which is the y-dependent effective mass governing the transverse momentum dependence. We obtain 


Uo(y,r T ) 

Ui(y,r T ) 


9 AM N ysJy j V M * K 0 {M T r T ) 

Mj 1 


2nFl 


where Kq and I\\ are the modified Bessel functions. At large values of the argument they behave as 




2 \J Mx'Tt 


{Mttt !)• 


(4.13) 

(4.14) 

(4.15) 

(4.16) 


The coordinate-space wave functions fall off exponentially at large transverse distances rx, with a width that is given 
by the transverse mass Eq. (14.141) and depends on the pion momentum fraction y. This behavior can directly be 
traced to the singularity of the momentum-space wave function at zero invariant mass, AM 2 = 0, which occurs at 
complex values of the transverse momentum, kj, = — Mjl, cf. Eq. (14.131) . 

The parametric domain in which we are interested in the coordinate-space wave function is 

y = 0(M n /M N ), r T = (4.17) 


In momentum space this corresponds to the region where the pion’s light front momentum components in the nucleon 
rest frame are 


k+ = yM N = 0(M n ), \k T \ = 0(M„), (4.18) 

and also k~ = (I^tI 2 + M 2 )/k + = 0(M n ), such that all components of the pion’s 4-momentum are 0(M V ) (“soft 
pion”). In this region chiral dynamics is applicable, and the approximations made in evaluating the peripheral 
contributions to the current matrix element are self-consistent. Equation (14.141) shows that for momentum fractions 
y = 0(MJM n ) 


M T {y) = 0(M n ) [y = 0(M n /M N )}, (4.19) 

so that the exponential range of the coordinate-space wave function is indeed of the order 0(M“ 1 ), cf. Eq. (14.161) . 

We note that for y = 0(1) the effective mass Eq. (14.141) is Mt = O(Mn), so that the range of the wave function 
Eq. (14.151) is 0(M^ 1 ). While the wave is still formally defined by Eq. (14. 1 5[) . it does not correspond to a chiral 
long-distance contribution in this case. This region does not contribute to the peripheral transverse densities, as the 
wave functions for y = 0(1) are exponentially small if the distance is kept at values rx = 0{M~ 1 ). In the calculations 
in Sec. |V]we can thus formally integrate up to y = 1 without violating the parametric restriction Eq. (14.171) . 

It is interesting to compare the parametric order of the light-front helicity-nonflip ( L z = 0) and flip ( L z = 1) 
components of the coordinate-space wave function in M^/Mn- Inspection of Eq. (14.151) shows that for y = 0{M n /Mj^) 
and r T = 0{M *.), Eq. (14.1711 , 


Uo/lh = 0 ( 1 ). 


( 4 . 20 ) 
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(a) ' (b) 

FIG. 4. Peripheral chiral light-front wave function in coordinate representation, (a) Radial wave functions Uq and Lh, Eq. (14.1511 
as functions of the pion momentum fraction y, at fixed transverse separation rr = l.OMj 1 . (b) Radial wave functions U\ as 
function of y at several transverse separations, tt = (1.0, 1.5, 2.0) M^ 1 . 


The helicity-nonflip and flip components are thus of the same order in the region of interest^ Regarding the numerical 
values we note that 

Uo(y,r T ) < Ui(y,r T ) (0 < y < 1, r T > 0 ), (4.21) 

because Mr(y) > yMu, cf. Eq. (14.1411 . and K\(z) > Kq(z) for all z > 0. The radial wave functions thus obey a 

numerical inequality at all values of the argument. 

Figure^] shows a plot of the peripheral radial wave functions f/o,i(y, rx) as functions of the pion momentum fraction 
y. Plot (a) compares Uq and U\ at a fixed transverse separation. One sees that the L z = 0 and L z = 1 components 
become equal at y —» 1 (i.e., at values several times M n /M N ), but show different power-like behavior at y — > 0, as 
is already apparent from the analytic formulas Eq. (14.1511 . One also sees that the inequality Eq. (14.211) is satisfied. 
Plot (b) shows the L z = 1 wave function U\ at several transverse separations. One sees that values of y ~ 1 are 

strongly suppressed with increasing transverse separation, and that the maximum of the wave function in y shifts 

to smaller values, in accordance with general expectations. At rr = several times M~ l the wave function and 
in particular the probabilities - are strongly concentrated at pion momentum fractions y = and the 

parametric approximations are borne out by the numerical results. 


V. PERIPHERAL TRANSVERSE DENSITIES 


A. Overlap representation 


We now want to express the peripheral transverse densities in the nucleon in terms of the chiral light-front wave 
functions. For this we first need to obtain explicit expressions for the invariant form factors in terms of the spin 
components of the current matrix element Eq. (ED. Taking the nucleon spin states in Eq. EH as light-front helicity 
states, cf. Eq. EH, and choosing a frame where the momentum transfer has only transverse components, cf. Eq. EH, 
the matrix element of the plus component of the current has the form 

(iV 2 | J + (0) \Ni) = (iV(p + ,p 2 T,cr 2 )| J + (0) \N(p + ,p 1T ,<Ji)} 


= (2 P + ) 


S(o -2,0-1) Fi(-A|) 


i (A t x e z ) ■ 5t(ct 2 ,cti) 


Mn 


(5.1) 


3 At exceptionally small pion momentum fractions y <C M^/Mn that the helicity-nonflip component of the wave function vanishes faster 
than the helicity-flip one, Uo/Ui —»• 0. This scenario is realized in the “molecular” region described in Ref. Q. 
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where St is defined as in Eq. (14.41) in terms of the rest-frame 2-spinors describing the initial and final nucleon, and 
e z is the unit vector in ^-direction. The form factors iq and iq are then obtained from the diagonal and off-diagonal 
matrix elements at0 


Fi(-A |0 ' 


2 ~+ \J + (0)\ N i) 

^ &1CT2 


-£(>1,(72) 


2 M n . , . 

-TT" H)( At X e ^) ‘ ^T(cri,CT2) 


(5.2) 


For the intermediate-nucleon part of the current matrix element we now substitute the overlap representation in 
terms of the light-front wave functions in the transverse rest frame, Eq. (14.101) . Using the coordinate representation 
the overlap integral becomes diagonal in the transverse relative coordinate rx and takes the form 

(N 2 \J + (0)\N 1 )y nteim = ^-^ J f d 2 r T e- lyrTA ‘ T ^2$*{y, r T ;a,a 2 ) <f>(y, r T;v,cr l)- (5.3) 


Notice that the momentum transfer A t is Fourier-conjugate not to rx itself but to yrx, which is a general feature of 
light-front kinematics. It is now straightforward to evaluate the spin sums in Eq. (15.21) and obtain the invariant form 
factors in terms of the L z = 0 and 1 components of the coordinate-space wave function Eq. (14.121) . We immediately 
quote the results for the isovector transverse densities p\ and 'p\ , Eq. m-- 

Pi (b) 1 

P%( b ) J 

The form of Eq. (|5.4I) is explained by the spin structure of the transitions (see Fig. [3 Jd). The light-front wave 
function has a nucleon helicity-conserving (Uq) and a helicity-flipping component ( U \). The current matrix element 
with the same nucleon helicity in the initial and final state requires the combination of two helicity-conserving or two 
helicity-flipping wave functions (Uq or Uf), whereas the matrix element with different nucleon helicities in the initial 
and final state requires combination of one helicity-conserving and one helicity-flipping wave function {UqU\). 

The transverse charge density p\ (5) also receives a contribution from the effective contact term in the current 
matrix element, Eq. m - This contribution cannot be represented as an overlap of nN light-front wave functions 
and has to be added to Eq. |QD as a separate term. The exact form of this term and its interpretation are discussed 
in Sec. IV FI The numerical contribution of the contact term turns out to be very small at distances b = few M~ 1 , so 
that the entire p\ is to good approximation given by the wave function overlap Eq. (15.41) . We may therefore compare 
the properties of the densities p\ and pX on the basis of Eq. (15.41) (the contact term is absent in p\ ). 


i [ d y j Vy)] 2 + 

yy 3 { -2U 0 (y,b/y)lh(y,b/y) 


(5.4) 


B. Chiral order and inequality 

The overlap representation Eq. (15.41) reveals several interesting properties of the chiral component of the peripheral 
transverse densities. First, because the light-front helicity-conserving and -flipping wave functions appear in the 
same order of the chiral expansion (their coefficients involve the same power of M n /Mix), cf. Eq. (14.201) . we conclude 
that the peripheral densities p\ and p\ are of the same order in the chiral expansion, 

(%{b)/pX(P) = 0(1) [b = 0(M~ 1 )}. (5.5) 

While this circumstance was noted earlier in the dispersive approach f8], where it is encoded in the chiral order of 
the spectral functions of the form factors near threshold, the overlap representation exhibits it more directly and 
provides a more immediate physical explanation. Notice that the original transverse magnetization density pX is 
parametrically larger than pX by a power of because the spatial derivative in Eq. (12.71) “counts” as 0(M 7T ) 

in the region b = 0(M~ 1 ), and thus 

(%(b)/pY(b) = 0(M n /M„). (5.6) 

The parametric equality of p\ and ~pX ■> Eq. (ESI), allows for non-trivial dynamical relations between the two densities 
and provides additional motivation for working with pJf rather than . 


4 The identification of the different spin components can be done conveniently by writing both sides of Eq. (15.ID as bilinear forms in the 
nucleon two-spinors, stripping off the two-spinors, and treating the equation as a 2 x 2 matrix equation. The different components can 
then be projected out by taking appropriate traces of both sides. 
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FIG. 5. Integrands of the transverse densities p\ 
and — pX in the wave function overlap representa¬ 
tion Eq. at a distance b = 1.0M, 1 . The 

integrands are given in units of M~ 2 . 


0 0.2 0.4 0.6 0.8 1 

Second, we observe an inequality between the peripheral transverse densities, 

\pX(b)\ < \pY(b)\- (5.7) 

It follows from the inequality obeyed by the quadratic forms in the integrands of the y-integral in Eq. (El, u'i + iff > 
2UqU\. While Eq. (15.71) was observed accidentally in the numerical calculations of Ref. | 8 j, its mathematical proof 
and physical explanation become possible with the wave function overlap representation provided here. Recalling 
that p\ ( b ) and pX (b ) represent the spin independent and -dependent components of the “plus” current at transverse 
position b in a nucleon localized at the origin and polarized along the y-direction, cf. Eq. (12.61) and Fig. [H we see that 
the inequality Eq. implies the positivity condition 

(J + (b))localized > \pX (b)\ - \pX (b)\ > o. (5.8) 

This property appears natural when one realizes that the plus current at b = 0(A1~ 1 ) results from peripheral pions, 
and that the current carried by an on-shell pion is proportional to its 4-momentum, ( 7 r + (fc)| J + | 7 r + (fc)) = 2 k + > 0. 
Such a “quantum-mechanical” picture of the peripheral densities will be explored further in Sec. IV El 



C. Numerical evaluation 

We now want to use the overlap representation Eq. (l5~fl) to study the numerical behavior of the transverse densities. 
Figure [5] shows the integrands of p\{b) and —p2(b) as functions of the pion momentum fraction y at a fixed transverse 
distance b = 1.0 M” 1 . One sees that the integrands are concentrated around values y ~ M n /M m- Contributions from 
y > 0.5 are very strongly suppressed because the wave function is evaluated at separations tt = b/(\ — y) that are 
substantially larger than b , and the wave function decays exponentially at large rr with a range that itself decreases 
with increasing y. One also sees that the integrands for the densities p\ and — p\ are close to each other at large values 
of y and differ only at y —> 0 , such that they are numerically close throughout the dominant region of integration. 
This follows from the similarity of the wave functions U$ and Ui, cf. Fig. HJi, and implies that p\ (b) « — /)\ (b) at 
distances b = few times M” 1 . 

The transverse densities obtained by performing the y-integral in Eq. (15.41) are shown in Fig. [ 6 ] One sees that 
/)X < 0 , and that the absolute value of the spin-dependent density is smaller than the spin independent one, | pX \ < p\ -, 
as required by Eq. UEZD- The inequality is almost saturated at distances b ~ 1 M n - 1 , as suggested by the integrands 
shown in Figure [5l but at larger distances \pX\ becomes significantly smaller than p\ ■ We note that the numerical 
densities obtained from the wave function overlap representation Eq. (15.41) exactly reproduce those calculated in the 
dispersive approach of Ref. [§}, which provides a test of the calculational procedures!! 


5 In Fig.8 of Ref. Q the function represented by the dashed line is — p 2 (b ), not p r z(b) (the plot is labeled incorrectly). The dispersive 
calculation gives p 2 (b) < 0, as does the wave function representation of the present work. 
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FIG. 6. LO chiral component of the nucleon’s 
isovector spin-independent current density p\ (6) 
(solid line) and spin-dependent current density 
—pX{b) (dashed line). The plot shows the densi¬ 
ties with the exponential factor exp(—2 M n b) ex¬ 
tracted [the functions plotted correspond to the 
pre-exponential factor in Eq. 


D. Transverse polarization 


Further insight into the peripheral transverse densities can be gained by considering the case of transversely polar¬ 
ized nucleon states. Transverse polarization naturally explains the similarity of the spin -independent and -dependent 
densities, p\ ( b) and pX (6), at distances b ~ few M/r 1 (cf. Fig. [6]) and enables a simple quantum-mechanical interpre¬ 
tation in the nucleon rest frame. 

Transversely polarized nucleon states in the light-front formulation are obtained by preparing a transversely po¬ 
larized state (say, in the p-direction) in the rest frame and performing the longitudinal and transverse boosts to the 
desired light-front momentum (cf. Sec. IIV All . With the general formula Eq. (14.11) this is accomplished simply by 
choosing the rest-frame 2-spinors as eigenspinors of the y-spin operator S v = cr y /2, 

Xtr (r = 1/2) = ^ ( • ) , Xtr(r = -1/2) = -±= ( l ) (5.9) 

(here and in the following we use r = ±1/2 to denote the y- spin eigenvalues). The vertex function calculated with 
the transversely polarized 4-spinors obtained in this way is of the same form as Eq. (Oil : only the structures S z and 
St = ( S x , S v ) are replaced by the contraction of the spin operators with the transversely polarized spinors 

= Xtr( T )(! CT *)Xtr(n) (i = x,y,z). (5.10) 

Notice that now the p-component is diagonal, 

>S? r (T,Tl) = t S(t, ti), (5.11) 


while the z and x components have off-diagonal terms. We define the light-front wave function for the N —> irN 
transition for nucleon states characterized by their transverse spin projections t\ and r, in complete analogy to 
Eqs. 02D gSD > 


^ti(y,k T -,T,T i) 


r t r(y,fcT;T,Ti) 


r t r(2/,fcr;T,n) 

AM 2 (y 1 k T ) 


2tg A M N 

FnVy 


yM N S z t (t , ri) ± k T ■ Str,T(r, n) 


(5.12) 

(5.13) 


Similar expressions are obtained for the nN —> N with transverse spin projections r and T 2 - 

The coordinate-space wave functions for transverse nucleon polarization are introduced through Eq. (14.111) in the 
same way as for longitudinal polarization and denoted by $ tr . The general decomposition Eq. (14.121) applies to the 
Fourier transform of the transversely polarized wave function as well, as it relies only on the functional dependence 
of the momentum-space wave function on kx, not on the specific form of the spin structures. Using the algebraic 
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relation between the spin structures for y— and z-polarization it is straightforward to express the transversely polarized 
coordinate-space wave function in terms of the longitudinally polarized radial wave functions Uq and U\. 


$tr(2/, r T , t = +1/2, n = +1/2) = sin aU\, (5.14) 

$tr (y, r T , t = -1/2, n = -1/2) = — sin aU\, (5.15) 

$tr(y, r T , t = +1/2, T\ = -1/2) = U Q + cos aU\, (5.16) 

$tr(y, r T , r = -1/2, ri = +1/2) = -U 0 + cosaUi, (5.17) 


r T = (rr cos a, 7 ’t sin a), t/ 0 ,i = Uo,i(y,rr)- 

With the radial wave function f7 0 and U\ given by the explicit expressions of Eq. (14.151) . these relations completely 
determine the transversely polarized coordinate-space wave function. 

In the transversely polarized representation the rotational symmetry around the z-axis is encoded in relations 
between the transverse spin components of the wave function. In the explicit formulas Eqs. (15.141) (15.171) these relations 

manifest themselves in that the four transverse spin components are expressed in terms of only two independent radial 

functions. One sees that the components satisfy 

$tr (y, r T , +, +) = -$tr (y, r T , -), (5.18) 

$tr (y, r T , +, -) = -$tr (y, -r T , +)• (5.19) 

The relation Eq. (15.191) between the transverse spin-flip components has a simple physical interpretation in the 
nucleon rest frame. In the wave function with initial nucleon transverse spin t\ = +1/2 and intermediate nucleon spin 
r = —1/2, the pion in the intermediate state has orbital angular momentum L = 1 with projection L y = +1 on the 
y— axis. Likewise, in the wave function with t± = —1/2 and r = +1/2, the pion has L v = —1. The two components 
thus differ only in that the pion rotates in the opposite sense around the y-axis, and one can be turned into the other 
by inverting the direction of the T-axis, i.e., replacing cos a —> — cos a in Eqs. (15.161) and (15.171) . 

The connection between the transversely and longitudinally polarized light-front wave functions generally depends 
on the angle of the transverse coordinate vector vt, as required by rotational invariance. A particularly simple 
connection is obtained at points on the negative or positive x-axis, where sin a = 0 and cos a = — 1 or +1. Considering 
the transverse spin-flip wave function on the negative and positive x-axis (“left” and ’’right” when looking at the 
nucleon in the z-direction from +oo, see Fig. |T|, and introducing the short-hand notation 

U\ e ft(y, r T ) = $tr (y, ~rre x , r = -1/2, n = +1/2), (5.20) 

bright(?/> r T ) = 4> tr (y, +r T e x , r = -1/2, n = +1/2), (5.21) 

we obtain from Eq. (15.171) 

Uieh(y, r T ) = ~U 0 (y, r T ) - Ux(y, r T ), (5.22) 

Uri g ht{y, r T ) = - u 0 (y , r T ) + U^y, r T ). (5.23) 

Note that the functions C/ieft aud I/ r i g ht are derived from a single y-spin component of the wave function but refer to 
a specific spatial direction. 

We can now derive a simple representation of the transverse densities in terms of the transversely polarized light- 
front wave function. Using the overlap representation Eq. m in terms of the the longitudinally polarized wave 
functions, and substituting them by the “left” and “right” transverse spin-flip wave functions according to Eqs. (15.221) 
and (15.231) . we obtair@ 

P\(b) | 

P%( b ) J 

This result can be explained easily by referring to the interpretation of the transverse densities as current matrix 
elements in a nucleon state localized in transverse space, cf. Sec. HT1 Fig. [ 1 ] and Ref. 0. According to Eq. (12.61) in 


1 r dy f [Uieft(y, h/y)] 2 + [f/right(y, b/y)} 2 | 
4jr VV \ -[{/ieft(y, b/y)} 2 + [{/right(y, b/y)} 2 j 


For simplicity we derive the representation Eq. (|5.24|) from Eq. (|5.4|) . using the relation between the longitudinally and transversely 
polarized wave functions at the special points b = ±be x - Equation (15.241) could equivalently be derived by converting the original current 
matrix element Eq. d5.ll) to the transverse spin representation and repeating the steps of Sec. IV Al In the latter approach one could 
choose any orientation of the vector 6; the rotational invariance of the densities would be guaranteed by the conditions Eqs. (I5.18D and 
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FIG. 7. Integrands of the transverse densities p\ 
and pX in the transversely polarized wave func¬ 
tion overlap representation Eq. (15.2411 . at a distance 
b = 1.0 M~ x . Solid line: [U\ e ft{y ,b / y )] 2 dashed 
line: [I/right (2/ , b/y)] 2 . The functions are plotted in¬ 
cluding the factor l/(4ivyy 3 ) and are given in units 
of M~ 2 . 
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a nucleon state with ?/-spin projection +1/2 (i.e., t\ = T 2 = +1/2) the dependence of the matrix element on the 
coordinate b is of the form 

(^ + ( b ))localized = (-) [Pl( b ) + COS 0^(6)], (5.25) 

where <f> is the angle of b relative to the i-axis. The densities are functions of b = |6|, and the angular dependence 
is given entirely by cos<f>. By choosing the direction of b along the positive and negative x-axis (cos </> = ±1), the 
individual densities can be expressed as 


(J + (be x ) + J+(-be x ))? oc aHzed = (■■•) pY(b), 
(J + (be x ) - J + (-be x ))Y ocalized = (...) p%(b). 


(5.26) 

(5.27) 


This structure is exactly analogous to Eq. (15.241) . with the integral over \U\ e {t(y, b/y)| 2 and |b/i e ft(y, b/y)| 2 representing 
the current densities at b = ±be x , respectively. 

The explicit form of the “left” and ’’right” transverse spin flip wave function, Eqs. (15.201) and (15.211) . is readily 
obtained from Eqs. (15.221) . (15.231) and (14.151) : 


U\ e h(y,r T ) 

bright iy i t t) 


9aM n yypy 

27 rF n 


\yMft Ko(MtTt) + MtKi(MtTt)] 

\yMj\[ Ao(H^tT't) — A/7 1 K\ ( AIt f'T )\ 


9 aM n Vy /y ( ( V M n + M t ) ) e -M T r T 
2V2\ ( y M N - M t ) J VMtTt 


(ttMt 1 ). 


(5.28) 


(5.29) 


The last expression is obtained with the asymptotic form of the modified Bessel functions, Eq. (14.161) . and applies at 
Mttt 1. Parametrically the two functions are of the same order in M^/Mys in the region y = 0 (M v /Mn) and 
rx = 0(M~ 1 ). Numerically one observes that 


\U\eit{y,r T )\ > |E/ri g ht(2/,rr)| {y ~ few times M V /M N ), (5.30) 

because Mriy) ~ yM n in this region of y. This is illustrated by Fig. [7] which shows the contributions of |C/i e ft| 2 and 
|bright| 2 to the integrands of the transverse densities in Eq. (15.241) . The strong suppression of the “right” compared 
to the “left” wave function is the reason for the similarity of the densities p\ ( b) and ( b ) at b = 0(M~ 1 ). 


E. Quantum—mechanical picture 

We can interpret our findings in a simple quantum-mechanical picture of peripheral transverse nucleon structure in 
the rest frame. For this purpose we imagine that the N -A nN transition takes place in ordinary time, and that the 
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initial/final state ^ •* intermediate state 


FIG. 8. Quantum-mechanical picture of peripheral transverse densities in LO chiral EFT. The bare nucleon with //-spin 
projection n = +1/2 in the rest frame (left) makes a transition to a pion-nucleon state with r = —1/2 and L y = +1 (right), 
and back to a bare state with T 2 = +1/2. The peripheral transverse densities p\ ( b ) and p%(b) are the left-right average and 
left-right asymmetry of the plus current (viewed from 2 = +oo) at the positions y = 0 and x = +6. 


wave function has the usual 3-dimensional rotational symmetry. For a non-relativistic system there would be a direct 
correspondence between the equal-time and the light front wave functions; see e.g. Ref. [24|. The chiral nN system is 
essentially relativistic, k = 0(M n ), and one should not expect a similar connection between the wave functions here. 
Nevertheless the intuitive quantum-mechanical picture of the chiral process explains all the essential features of the 
peripheral transverse densities. In any case its content is backed up by light-front wave function formulas, which are 
exact also in the relativistic case. 

Consider a nucleon in the rest frame, in a spin state polarized in the positive (/-direction (S v = +1/2). In the 
interaction picture implied by chiral EFT, we may think of this physical nucleon as a pointlike bare nucleon coupled 
to soft pions, described by a wave function. The LO contribution to the peripheral charge and current densities in 
the proton isospin state arises from the component with a single peripheral positively charged pion. This component 
corresponds to the chiral process where the initial (bare) proton makes a transition to a state with a (bare) neutron and 
a peripheral positive pion, and back to the final bare nucleon (see Fig [5]). The y- spin projection quantum numbers of 
the initial/final nucleon state are t\ = T 2 = +1/2, and that of the intermediate nucleon state is denoted by r. Because 
of parity conservation the wave function of the 7 r + n system has orbital angular momentum L = 1. Conservation of the 
total angular momentum projection on the y— axis allows only the states with L v = 0 and r = +1/2 (spin-conserving), 
and with L v = +1 and r = —1/2 (spin-flip). We are interested in the densities in the x—z plane (y = 0), to which 
the state with L v = 0 cannot contribute, as its wave function vanishes in the direction perpendicular to the //-axis, 
Pi (cos 9) = 0 for 9 = 7 r /2 (Pi is the Legendre polynomial of degree 1). This leaves the state with L v = +1 and 
r = —1/2 as the only contribution to the densities in question It explains why in the light-front formulation we were 
able to express the transverse densities completely in terms of the transverse spin- flip wave function, cf. Eqs. (15.201) 
and (15.211) . 

The peripheral densities in the x—z plane thus arise from configurations in which the positive pion “orbits” around 
the neutron with angular momentum L v = + 1 . Since the pion itself has no spin, the current it produces is the 
convection current caused by the orbital motion of the charge. Because the chiral EFT interactions between the 
pion and the nucleon have a short range -C M ~ 1 , the peripheral pion can be regarded as free while the current 
is measured. The 4-vector current carried by a free pion with momentum k and charge +1 is = 2 k^, where 
k° = Ett = \/|fc | 2 + M 2 , so that its plus component is positive for all momenta, J + = 2 k + = 2 (E n + k z ) > 0. This 
explains the positivity property of the current density, Eq. m - 

Now according to Sec. |TT] and Eq. (12.61) the transverse densities p\ and pi/ are the left-right average and left-right 
asymmetry of the J + current density produced by the peripheral pion. It is obvious that a positively charged pion 
orbiting with L v = +1 produces a larger J + density on the left (where it moves in the positive z direction) than on 
the right (where it moves in the negative z direction). This explains why p \(b) < 0 (see Fig. [6]). The magnitude of 
the asymmetry is determined by the effective pion momenta and the relativistic effects implied by the projection on 
fixed light-front time. If the motion of the pion were non-relativistic, with characteristic velocity v = k/M w -C 1, the 
plus momentum carried by the pion would be k + = M n [l + 0(i>)]; i.e., it would be dominated by the pion mass and 
independent of the direction of the pion momentum. Since furthermore the probability to find a pion would be the 
same on the left and on the right side of the x axis (because of rotational symmetry around the y-axis) the ratio of 
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left and right current densities at the same distance b = \b\ would be 


(T^( localized 

localized 


1 + 0(v). 


(5.31) 


The light-front wave function representation of Sec. IV Dl shows that the asymmetry is much larger than the non- 
relativistic estimate Eq. (I5.31|) . 


(T~*"( 5ea;))localized 

( V^da;)) localized 


I ^ risht ^ y ’ ^)l 


(5.32) 


The numerical value of the ratio is ~ 9 at b = 1 M~ x and ~ 4 at b = 5 M~ x (see Fig. [TJ). This highlights the 
essentially relativistic nature of the motion of pions in chiral dynamics. The power of the light-front formulation is 
that it permits a first-quantized representation even of such essentially relativistic systems. 


F. Contact term 

The transverse charge density p\' (b) also receives a contribution from the effective contact term in the current 
matrix element, Eq. m- This contribution cannot be represented as an overlap of nN light-front wave functions 
and needs to be added separately to that from the intermediate nN state, Eq. (E3D, 

Pi (b) = Pi (fr)interm [from Eq. (|5.4|)] + p\ (6) contact • (5.33) 

It is readily computed by evaluating the Feynman integral Eq. (13.911 as a four-dimensional integral, using the fact 
that it depends only on the momentum transfer A M as external 4-vector. One obtains Q 

pX(b) contact = (1 192 $^ {[K 2 (M„b)] 2 - 4[K 1 (M n b)] 2 + 3 [K 0 (M„b)] 2 } . (5.34) 

The contact term density is negative because 1 —g\ < 0. It is of the same parametric order as the one from intermediate 
■nN states, as can be seen by comparing Eq. (15.3411 with Eqs. (15.411 and (14.1511 and noting that at distances b = 0(M~ 1 ) 
the integral is dominated by pion momentum fractions y = 0(M 7r /M/v). The numerical contribution of the contact 
term density turns out to be very small in the region of interest, b = few M ranging from about —10% of the 
intermediate nN contributions at b = 1M” 1 to —4% at b = 5 M~ l (see Fig. [9]). Thus the entire p\ is to good 
approximation given by the wave function overlap Eq. (El, which justifies our earlier comparison of the properties 
of p\ and p\ on the basis of Eq. (15.411 (the contact term is absent in pX )■ 

Some comments are in order regarding the interpretation of the contact term in the context of the light-front 
description. In the intermediate nN contribution to the current matrix element the typical energy denominators are 
AM 2 (y,kT,piT, 2 T) = 0(M ot Mjv), cf. Eqs. (13.1811 and (13.1911 . and the typical pion light-front energies in the nucleon 
rest frame are fcj" 2 = 0{M 1r ). These are configurations that exist for a large light-front time interval Ax~ = 0(M ~ 1 ) 
and can be regarded as particle states in the chiral EFT. The contact term describes contributions to the peripheral 
density from intermediate states with invariant mass differences that are not chirally small, i.e., that do not vanish in 
the limit M n —> 0. These states lie “outside” the chiral EFT, and their contribution is represented by a local operator. 

The appearance of the combination 1 — g\ in the coefficient of the local operator is natural |8(. For a nucleon 
without internal structure one would have gA = 1, and the effective contact term would be absent. The deviation of 
gA from unity is the result of the “compositeness” of the nucleon, which in turn is related to the presence of inelastic 
states in nN scattering. The combination 1 — g\ can thus be regarded as the effect of non-chiral intermediate states 
in the nN scattering amplitude, in agreement with the above interpretation^ Further properties of the contact term, 
such as its relation to the form of the nNN coupling (pseudoscalar, axial vector) are discussed in Ref. Q- 

The contact term formally corresponds to a contribution of light-front zero modes (vanishing pion plus momentum) 
to the current matrix element. An advantage of our approach, starting from invariant integrals, is that it allows us 
to identify and calculate these contributions in a straightforward manner. The zero mode contributions could also 
be calculated in the time-ordered formulation of chiral EFT, by considering the current matrix in a frame where the 
current transfers plus momentum, A + ^ 0, and taking the limit A + —> 0 at the end of the calculation. 


' The connection between g\ — 1 and inelastic states in nN scattering is expressed in general terms by the Adler-Weisberger current 
algebra sum rule l2ol . 12(4 . 
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b [M* 1 ] 


FIG. 9. LO chiral component of the nucleon’s 
isovector spin-independent current density pXib). 
Solid line: Contribution from intermediate nN 
states, Eqs. El and (lOll . Dotted line: Effec¬ 
tive contact term, Eq. (13791) . plotted with oppo¬ 
site sign (the actual contribution is negative). The 
plot shows the contributions to the density with 
the exponential factor exp(—2A/ 7r f>) extracted [the 
functions plotted correspond to the pre-exponential 
factor in Eq. El]- 


VI. CHIRAL GENERALIZED PARTON DISTRIBUTIONS 


A. Peripheral pion distribution 


In QCD the transverse structure of the nucleon is expressed in terms of coordinate-dependent parton densities, 
which are defined as the Fourier transforms of the generalized parton distributions and describe the density of partons 
(quarks, antiquarks, gluons) with a given light-front plus momentum fraction x at transverse position b. In this 
context the transverse charge and current densities are obtained as integrals of the transverse densities of charged 
partons (quarks minus antiquarks) over x. Chiral dynamics governs not only the peripheral charge and current 
densities but also the x-dependent distributions of partons at b = |6| = 0(M~ 1 ) and x = O(/ M,.y). Detailed 
studies of the peripheral parton densities due to chiral dynamics have been performed in Refs. [13, [l8|. Here we 
show that the x-integral of these peripheral parton densities reproduces the transverse charge and current densities 
calculated in chiral EFT. We also express the plus momentum distribution of peripheral pions in the nucleon in terms 
of the nN light-front wave functions introduced in Sec. IIIII 

The basic object in the study of peripheral partonic structure is the light-front plus momentum distribution of soft 
pions in the nucleon. Following Refs. [313 we define the GPDs of soft pions in the nucleon in terms of the matrix 
element of the bilinear light-ray operator in the pion field 


d£_ 

2n 


iyp + £ / 2 


{N(p 2 ,a 2 )\ 


£ e 3ab n a (-£/2) d + n b (£/2) | €+=(Ur=0 


_ ab 


|N(pi,cti)) 


= U 2 


+ ZJV, ^ E n(y^) 

7 Hit (V, t) -~- er 


2Mn 


u 1 , 


( 6 . 1 ) 


where £ is the 4-vector of the space-time separation of the fields and = (<9^— <9 M )/2. Equation (16.11) applies 
in the parametric regime of pion momentum fractions y = 0 (M v /Mn) and momentum transfers t = 0(M%) and is 
to be evaluated in chiral EFT. The pionic operator on the left-hand side has isovector quantum numbers, and the 
matrix element is understood in the sense of Eq. (12.81) . The GPDs are denoted by H\ and in accordance with 
the standard convention [23], and are defined in the interval — 1 < y < 1. It is easy to see that the pionic light-ray 
operator is symmetric under £ —> — £, whence the isovector GPDs are even functions of y, 


= H^(-y,t), 


E^{y,t) = E^{-y,t). 


( 6 . 2 ) 


Note that for £ = 0 the pionic operator reduces to the local vector current of the pion field, 

J+(0) = e 3ab 7T a (0), a+7T & (0), 

ab 


(6.3) 
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the matrix element of which determines the chiral 7 T 7 t cut contribution to the Dirac and Pauli form factors and was 
calculated in Sec. IIII Al When integrating Eq. (16.11) over y the exponential factor produces a delta function which 
enforces = 0 and thus £ = 0 , so that the matrix element of the light-ray operator becomes that of the local vector 
current. In this sense Eq. ED is just a particular representation, differential in the pion plus momentum fraction y , 
of the chiral 7 T 7 t cut in the current matrix element. 

The transverse coordinate-dependent (or impact parameter-dependent) distributions of soft pions in the nucleon 
are then defined in analogy to the transverse charge and current densities, Eq. (USD , as 


HX(y,t = - A*) 
EX(y,t = - A?,) 


/ 


d 2 b e iArb 


f¥*(y,b) ) 


(6.4) 


It is convenient to introduce a modified helicity-flip distribution in analogy to p 2 , Eq. (ED, 



2 M]y 


(6.5) 


The functions b) and /^(y, b) describe the isovector transverse spatial distribution of pions with plus momentum 

fraction y and refer to the parametric regime b = and y = The interpretation of these spatial 

distributions in a transversely polarized nucleon is analogous to that of the transverse densities pY (b) and pX(b), 
cf. Eq. (003) and Fig. □ 

The peripheral pion distributions /^(t/, 6 ) and /^.(y, b) can be calculated in LO chiral EFT in the same manner 
as the transverse densities p\ ( b ) and p\ ( b ), cf. Sec. IIIII One expresses the matrix element in Eq. (16.11) as a Feynman 
integral and separates it into the intermediate nucleon contribution and an effective contact term, cf. Eqs. ED 
ED- The intermediate nucleon contribution can then be represented as the product of the light-front wave functions 
describing the transition of the initial and final nucleon to the nN intermediate state. For y > 0, 

fin(V’ b ) 1 1 j [Uo{y,b/y)] 2 + [Ui(y, b/y)] 2 \ 

fl(y,b) J 2ny y 3 \ -2U 0 (y,b/y) Ui{y,b/v) J 

where t/ 0 ,i are the coordinate-space wave functions defined in Eqs. (14.121) and 
[cf. Eq. (16.21) ] 

fi*(y, b ) = fi*(-y> b )> fLiV’b) = r (-y,b). ( 6 - 7 ) 


(:y > 0 ), ( 6 . 6 ) 

(14. 15p : for y < 0 one uses that 


The contact term contribution is obtained from the Feynman integral as 

/&(», &)contact = S(y) pX {b) contact, (6.8) 

where pX(b) contact is given by Eq. (15.341) . This expression clearly identifies the contact term as a light-front “zero 
mode” contribution. The complete isovector pion distribution fY n is then given by the sum of the overlap contribution 
Eq. (16.61) and the contact term Eq. (16.81) . This representation of the chiral pion densities in particular implies that 
[cf. Eqs. (15.41) and (15.331) ] 



as is obvious from the above definition of fY and fY n - In this sense our earlier results for the peripheral transverse 
densities in chiral EFT have a straightforward interpretation as y-integrals of the peripheral pion GPDs in the 
nucleon. The significance of this connection lies in the fact that the peripheral pion GPDs have a more general 
physical significance and can in principle be measured independently in peripheral high-energy scattering processes. 


B. Charge density from peripheral partons 


We can now demonstrate the connection of the chiral component of the peripheral transverse densities with the 
peripheral quark/antiquark content of the nucleon in QCD. Following Ref. [181] the isovector quark/antiquark density 
in the nucleon at b = 0(AI~ 1 ) generated by chiral dynamics is given by 


= [d — u] (x, b) 


chiral — 


dy f V , n 

— fiAy, b ) y 
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( 6 . 10 ) 
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where z = x/y represents the fraction of the pion plus momentum carried by the quark/antiquark, and q^(z) is the 
valence quark/antiquark density in the pion, 

= ± [d-u\ n± (z) = ±[u-d\ 1T± (z) = ±§ [u-u- d + d\ n± (z), (6.11) 

normalized such that 


l dzq™\z) = 1. (6.12) 

Equation (16.101) has the form of the usual partonic convolution formulas and relies on the approximation that the 
non-chiral transverse size of the pion can be neglected on the scale i.e., the spatial distribution of peripheral 

quarks/antiquarks is determined entirely by the distribution of pions in the nucleon. The transverse charge density 
in the proton is generally given by the integral of the quark minus antiquark densities over X, weighted by the quark 
charges (e u = 2/3, e d = -1/3), 


Pi(b) = / dx {e u [u - u](x,b) + e d [d-d](x, b)} . 
Jo 


(6.13) 


The isovector density is obtained by taking half the proton-neutron difference and using isospin symmetry 


Pi(b) = b[Pi ~ PiK b ) = b(e u -e d ) dx [u - u - d + d\(x,b) (e u - e d = 1) 


(6.14) 


(this expression is valid even when including strange quarks in the individual proton and neutron densities). Sub¬ 
stituting here the peripheral quark/antiquark densities generated by chiral dynamics, Eq. (16.101) . and using the 
normalization condition Eq. (16.121) . one obtains 


p\ (fr)chiral = dx [u - U - d + d\(x, 6) c hiral 


= dx 


C\x/y) 


[ dyf?„(y,b) [ —q™\x/y) 
Jo Jo y 


= / dyf{ n (y,b), 


(6.15) 


which agrees with Eq. (16.91) . It shows that resolving the peripheral pion into its quark/antiquark constituents and 
computing the charge density from the quark/antiquark densities in the nucleon leads to the same result as computing 
the charge density directly from the distribution of (pointlike) pions — as it should be in the parton picture. The 
same applies, of course, to the spin-dependent density p% ( 6 ). 

The x — > 0 limit in the integral Eq. (16.141) and the role of the contact term in the charge density in the “sum rule” 
Eq. (16.151) require special consideration. In the LO chiral expansion of the current matrix element and the correlator 
Eq. (16.11) it is supposed that the pion plus momentum fraction is of the parametric order y = 0(M n /Mff). The 
contributions from y —> 0 on this scale are described by a delta function S(y), i.e., one sees only their contribution 
to the total charge density ('(/-integral, first moment) but cannot resolve their dependence on y. In this sense the 
LO chiral expression for the peripheral parton densities is valid for x = 0(M w /Mjv) but otherwise not exceptionally 
small. Now calculating the charge density requires integration down to x —> 0. It is clear from the above that the 
LO chiral expansion and the limit x —> 0 do not commute. It explains why Eq. (16.151) captures only the non-contact 
(intermediate irN) contribution to the charge density, and why the contact term has to be added separately. The LO 
chiral expansion correctly accounts for the total charge density; it is just not smooth enough at y —> 0 to allow for 
the charge to be distributed over partons with finite x and recovered by integration over x. A resummation of chiral 
EFT at parametrically small x in the logarithmic approximation has been proposed in Refs. (28l-[30j; this approach 
“resolves” the delta function of the finite-order approximation into a finite-width function that allows integration 
over x. We emphasize that the small-cc problem discussed here is largely formal, and that chiral dynamics likely 
does not dominate the actual small-a; behavior of peripheral parton densities (for a discussion of chiral dynamics and 
resummation in the light of conventional small -x physics, see Ref. (l 8 |). 














24 


VII. SUMMARY AND OUTLOOK 

Using the light-front representation of relativistic dynamics we have expressed the LO chiral EFT results for 
the nucleon’s peripheral transverse densities in “first-quantized” form, as overlap integrals of the light-front wave 
functions of a chiral irN system. The new representation is exactly equivalent to the “second-quantized” field- 
theoretical results and enables an intuitive understanding of chiral dynamics in close analogy to non-relativistic 
quantum mechanics. It reveals an inequality between the spin independent and -dependent transverse densities, 
\f> 2 (&)| < p\ (b). which constrains the spatial distribution of charge and magnetization in the nucleon. It also offers 
a simple dynamical explanation why the inequality is almost saturated. The wave function representation permits 
straightforward numerical evaluation of the transverse densities and the underlying plus momentum distributions. 
It also connects the peripheral transverse charge density with the nucleon’s peripheral partonic content (GPDs) 
generated by chiral dynamics. 

Our studies reveal two interesting general aspects of chiral dynamics and peripheral nucleon structure. One is the 
essentially relativistic character of chiral dynamics in the parametric region of pion momenta k = O(M^). The large 
left-right asymmetry of the chiral densities in the transversely polarized nucleon is a genuine relativistic effect and 
results in a ratio very different from the non-relativistic estimate. It predicts the approximate equality |(6)| ss p\ (6) 
at b = few M“, which can be tested experimentally. In this sense measurements of the nucleon’s Dirac and Pauli 
form factors can directly attest to the relativistic nature of chiral dynamics in the nucleon’s periphery (for a discussion 
of the prospects for extracting the peripheral transverse densities from nucleon form factor data, see Refs. urn). 

The other interesting aspect is the role of the pion’s orbital angular momentum in peripheral nucleon structure. It 
is seen most clearly in the case of transverse nucleon polarization, where a single pion orbital with L = 1 accounts 
for both the spin-independent and -dependent densities and explains their properties. This places chiral dynamics in 
the context of contemporary studies of orbital angular momentum in relativistic systems and quantum field theory, 
inspired by the nucleon spin problem (for a recent review, see Ref. [3l| ). The first-quantized light-front representation 
is an essential tool in defining the angular momentum content of relativistic systems and interpreting the dynamics, 
and is therefore natural for chiral dynamics. 

The chiral two-pion exchange contribution studied in this work affects the isovector component of the nucleon 
charge densities. In the isoscalar component the chiral contribution starts with three-pion exchange and is strongly 
suppressed, because of its shorter range and its higher order in the chiral expansion. In this sense the isovector 
component computed here determines the large-distance behavior also of the individual proton and neutron densities, 
which are the sum of isovector and isoscalar components. The analysis of experimental data for the proton and 
neutron form factors at low |t| should be done with dispersion-based parametrizations [32, 33j, which have correct 
analytic properties (singularity structure) and smoothly combine the two-pion exchange contribution with the vector 
meson resonances determining the bulk of the transverse densities [Hj]. 

The methods developed in the present work can be applied to several related problems in nucleon structure. Of 
particular interest would be the study of chiral dynamics in the matrix elements of the energy-momentum tensor, 
whose transverse densities describe the spatial distributions of matter, momentum, and stress (or forces) in the 
nucleon. This would in particular allow one to confront the “particle-based” definition of orbital angular momentum 
in the light-front representation with the “field-theoretical” definition in terms of the energy-momentum tensor, and 
in this way test the various angular momentum sum rules proposed in the recent literature [3514391 ] . 

The light-front wave function representation of peripheral densities could also be extended to A isobar intermediate 
states. The A contribution to the transverse densities was computed in the dispersive approach in Refs. Id. While 
is is numerically small at distances b = few M~ Y it plays an important role in ensuring the proper scaling behavior of 
the peripheral densities in the large-iV c limit of QCD, where the N-A mass splitting scales as Ma — Mjv = 0(N~ 1 ). 
The light-front wave functions for the N —> n + A transition can be defined in analogy to the N —> n + N ones, 
cf. Eq. (13.2611 . Light-front time-ordered calculations with higher-spin particles are generally plagued by ultraviolet 
divergences resulting from the breaking of rotational invariance. An advantage of the Lorentz-invariant approach 
taken in the present work (see Sec. ED, in which the light-front representation is derived from the Feynman integrals, 
is that it maintains rotational invariance and avoids such divergences. It is thus particularly suited to including A 
intermediate states. 

The nucleon’s chiral component can be probed also in high-energy scattering processes in yiV, eN, irN , or NN 
scattering (squared center-of-mass energies s>l GeV 2 ), by selecting reaction channels and kinematic regions where 
the process happens predominantly on a pion at transverse distances b = 0(M~ 1 ). Under certain conditions the 
amplitude for such processes can be expressed in terms of the light-front wave functions of the peripheral nN system 
and the amplitude for the high-energy scattering process on the pion. These are the same light-front wave functions 
as those introduced in the present study of peripheral current matrix elements. Only the light-front formulation of 
chiral dynamics makes it possible to establish such a connection between low-energy and high-energy processes. This 
new connection greatly enlarges the number of experimental probes of chiral dynamics. 
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One example is hard exclusive electroproduction of mesons on a peripheral pion, e + N — > e 1 + 7r + meson + N' at 
Q 2 1 GeV 2 and x -C M^/Mn, in the region where the invariant momentum transfer between the initial and final 
nucleons is tNN' = (pn—p'n) 2 = 0(M 2 ). Here the high-energy process on the pion, e + 7r —> e'+meson+7r, probes the 
quark/gluon GPDs in the pion E3- Such measurements could be performed at a future Electron-Ion Collider (EIC) 
with appropriate forward detectors for the pion and the recoiling nucleon. Another example is wide-angle quasi-elastic 
scattering on a peripheral pion, n +N tt+ n +N', in the region where tNN 1 = 0(M 2 ); here the high-energy process 
is described by the n + n —> n + tt elastic scattering amplitude at squared center-of-mass energies s n7r 1 GeV 2 . 
An interesting consequence of chiral dynamics is that the cross sections for such “pion exchange” processes should 
have large transverse single-spin asymmetries, governed by the transverse spin dependence of the chiral nN light- 
front wave function. This circumstance might further help to distinguish such processes from conventional “vacuum 
exchange” processes. 
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Appendix A: Light—front time—ordered formulation 


In Sec. m the light-front wave function of the chiral nN system is introduced as an element of a particular 3- 
dimensional reduction of the Feynman integrals for the nucleon form factor in relativistically invariant chiral EFT. 
In this appendix we show that this object is identical to the conventional wave function, defined as a transition 
matrix element in light-front time-ordered perturbation theory El- The correspondence is useful for relating our 
results to more phenomenological applications of light-front time-ordered perturbation theory, and for an eventual 
time-ordered formulation of chiral EFT. 

In light-front quantization we consider the evolution of chiral effective field theory in light-front time x + = t + z. 
It is governed by the Hamiltonian H = P~/ 2, where P = P° — P z is the minus component of the energy- 
momentum 4-vector of the field theory. In the interaction picture the Hamiltonian is split into a free part and an 
interaction, H = Hq + H m t. The eigenstates of the free Hamiltonian (which we denote as |...)o) are nucleon and pion 
single-particle states, characterized by their light-front plus and transverse momenta, with eigenvalues given by the 
light-front energies, 


and are normalized as 


Ho |N(p,cr))o 
H 0 |7r(fc)) 0 


Tj" \H(p, cr ))oj 
K(*0)o, 


p 


k~ = 


= Pt + m n 

2 p+ 

k 2 T + Ml 
2 k+ ’ 


o(N(p\a')\N(p,a)}o = 2 p + (2 tt ) 3 5(p' + - p + ) 6 {2) {p' T - p T ) S(a',a), 
o(7r(fc , )|7r(fc)) 0 = 2 k + (2n) 3 S(k ,+ — k + ) 5^ 2 \k' T — fey). 


(Al) 

(A2) 

(A3) 

(A4) 


The interactions represent an operator in the product space spanned by the single-particle states and induce transitions 
between product states with different particle number. The evolution of a state at light-front time x+ to is described 
by the action of the time evolution operator 


S(x 2 ) = T exp 


—i 


dx + H int (x + ) 


H int (x + ) = exp (iH 0 x + ) H int exp (~iH 0 x + ), 


(A5) 

(A6) 


where the time ordering is in x + and it is assumed that the original interaction operator H m t does not explicitly 
depend on x + . 

Consider now the matrix element of the current operator at time x + = 0 between physical nucleon states with plus 
momenta pf = p% = p + and transverse momenta such that P2T — Pit = A t , cf. Eq. (13.101) . Assuming that the 
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interaction is switched off adiabatically at x + —> ±oo, the matrix element is given by 


(N 2 \ J+(0) |fVr) = (N(p 2 ,a 2 )\ J+(0) ^(pi,^)) 

= o(-^(P2,ct 2 )| 5(oo,0) J+(0) 5(0, -00) |-/V(pi,fTi)) 0 . 


(A7) 


A peripheral contribution originates from the components of the evolved nucleon state that contain a peripheral pion. 
In leading order of the interaction H lnt this component is a state with a nucleon and a single peripheral pion. The 
corresponding contribution to the matrix element is obtained by inserting free pion-nucleon intermediate states to 
the left and right of the current operator, 


f dkt 

f d 2 kir 

f dk+ 

j ' d 2 k 2 T 

f dl * , 

f d 2 lr 

1 (2tt)2 kt J 

(27 r ) 2 J 

(2tt)2 k+ J 

(27T) 2 J 

(2tt)2Z+ J 

(2t) 2 


x 0 <JV(P2, )| S (°°’ °) \n(k 2 )N(l, a)) 0 

x oMMI -7 + (0) k(fci)) 0 


x 0 (n(ki)N(l, a )| 5(0, -oo) |iV(pi, cti)) 0 . 


(A8) 


The amplitude for the transition from the initial free nucleon state to the pion- nucleon intermediate state at time 
x + = 0 is 


o{ir(ki)N(l, cr)| 5(0, -oo) LO |lV(pi, cti)) 0 
r° 

= o(^(ki)N(l,a)\ (-£) / dx + H int (x + )\N(p 1 ,ai)) 0 

J — OO 

= o(n{ki)N(l,*)\H int \N(p 1 ,tT 1 ))o (-*) f dx+ e *K-*)*+/* 

J — OO 


0 (7r(fci)jV(l, cr)\ H int | N(px, q 1 )) 0 
h( k i + l ~ - Pi) 


Pi = 


Pit + fc— _ k iT + 


’ 'M 
Pi 


K 


1 2 t + M 2 n 
1+ 


(A9) 


where the infinitesimal imaginary part of the light-front energies implies that the interaction is switched on adiabat¬ 
ically as x + increases from — oo. The interaction Hamiltonian is the integral of the interaction Hamiltonian density, 
which in turn is the negative of the interaction Lagrangian density, 


H nt = 


2 J dx d^Xx T^int(^r)x+—0 — 2 J d^- d 


:+= 0 - 


We use translational invariance to evaluate the matrix element in Eq. (1X91) . 


(A10) 


o (7r(A;i)iV(Z,0-)|-Hi n t |lV(pi, <ti)) 0 

= ~ J dx~d 2 x T o{Tr(ki)N(l,a)\ C int (x) |lV(pi, oi)) 0 

= ~ J dx~d 2 x T e l(kl+l ~ Pl)x 0 (n(ki)N(l, a)\ £ int {0) \N(p 1 ,ai)) 0 

= —(2 tt 3 ) S(k+ +1+- p+) S^(k 1T + l T - p 1T ) 0 (Aki)N(l, a)\ £ int (0) |iV( Pl , ar)} 0 . (All) 

The LO transition matrix element Eq. (1A91) can thus be represented as 

0<7T (h)N(l, a)\ 5(0, -oo) LO \N( Pl , a^o = 2 p+ (2tt 3 ) S(k+ + l + - pf) S^ 2 \k 1T + l T - p 1T ) T, (A12) 




T(y, fc T ,PiT;cr, (Ji) 


o (7r(fci)lV(Z,g)|£ in t(0) |lV(pi, o~i))o 
PiiK + l ~ ~Pi) 


(A13) 


where T is the light-front wave function of the nN component of the nucleon. It is invariant under longitudinal 
boosts and can be regarded as a function of the independent momentum variables y = k+jp\ , kx = kyp ~Pit = — It , 
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and Pit■ An analogous expression can be written for the complex conjugate matrix element describing the transition 
from the outgoing free pion-nucleon intermediate state to the final free nucleon state. 

The wave function defined by Eq. (IA13I) is identical to the one introduced in the context of the reduction of the 
Feynman integral, Eq. (13.261) . Namely, the light-front energy denominator in Eq. (IA13I) is just the invariant mass 
difference Eq. (13.171) . 

pt(K + l- - pr) - 1T — + T +_y N - Ml - p\ T = AM 2 . (A14) 

Furthermore, the transition matrix element of the Lagrangian density in the numerator in Eq. (IA13I) is just the LO 
7 tN vertex function Eq. (13.231) . 

o(ir(ki)N(l,<j)\C int (0) \N(p 1 ,a 1 )) 0 = F(y, fc T , pi T ; cr, cti) + terms oc (fcf + l~- _pj“). (A15) 

The expression quoted here is for the p —> n 0 + p isospin component of the matrix element; the other components 
follow from isospin invariance [cf. the comments after Eq. (13.271) ]. 

The final result for the leading-order peripheral contribution to the current matrix in the time-ordered approach 
is then obtained by (a) substituting the transition matrix elements in Eq. (IASI) by the representation Eq. (IA12I) , with 
the isospin factor y/2 for the transitions p —> 7r + + n and n —> + p; (b) substituting the explicit expression for the 
matrix element of the current operator between free (pointlike) charged pion states, o( 7r ^ ±) (^ 2 )| J + (0) 17r (±) (*i))o = 
±(fcj l 11 " + Af]~); (c) integrating over the redundant intermediate pion and nucleon momenta using the delta functions 
in Eq. (IA12I) . The result for the isovector component is identical to Eq. (13.271) . Altogether this establishes the 
correspondence between the time-ordered and the invariant calculation of the LO chiral contribution to the peripheral 
densities. 

Note that the transition matrix element Eq. (IA15I) calculated with the axial vector ttNN coupling (as appears in 
the original chiral Lagrangian) and the pseudoscalar coupling (as emerges from our reduction of the Feynman integral) 
differ by a term oc (fcj~ +l~ —Pi). This term cancels the energy denominator of the wave function and therefore results 
in a contact term in the matrix element, in agreement with the findings of Sec. IIIII One thus understands how the 
effective contact term Eq. (13.91) appears in the time-ordered formulation. Calculation of the complete contact term 
contribution in the time-ordered formulation would be possible with a careful limiting procedure for the instantaneous 
exchanges. An advantage of our invariant formulation is that it allows us to calculate this contribution with minimum 
effort. 
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